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ABSTRACT 

We present new optical and infrared (IR) observations of Cir X-1 taken near apastron. 
Both sets of spectra show asymmetric emission lines. Archival optical observations 
show that an asymmetric Ha emission line has been in evidence for the past twenty 
years, although the shape of the line has changed significantly. We present an eccentric 
(e ~ 0.7-0.9) low mass binary model, where the system consists of a neutron star 
orbiting around a (sub-)giant companion star of 3-5 Mq. We suggest that the broad 
components of the emission lines arise in a high-velocity, optically thick flow near 
the neutron star; while the narrow components of the optical and the IR lines arise 
near the companion star and a heated ejecta-shell surrounding the binary respectively. 
In this model, the velocity of the narrow component reflects the space velocity of the 
binary; the implied radial velocity (-1-430 kms~^ after correcting for Galactic rotation) 
is the highest velocity known for an X-ray binary. 
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1 INTRODUCTION 

Cir X-1 is one of the most puzzling X-ray binaries known. 
Like the peculiar systems SS 433 and Cyg X-3, it cannot 
easily be classified into any of the major categories of X-ray 
binaries. Indeed, there is even doubt as to whether it is a 
high- mass (HMXB) or low- mass X-ray binary (LMXB). 

Since its discovery in early 1970s, Cir X-1 has been stud- 
ied intensively at X-ray wavelengths. The X-ray properties 
of Cir Jf 



from the X-ray binary. Fender et al. (1998) have imaged an 



arcsec-scale asymmetric jet aligned with these larger struc- 
tures, raising the possibility that the outflow from th e sys - 
tem is relativistic. Recently, Case & Bhattacharya (199J) 
have revised the estimated distance to G321.9— 0.3 (and 
hence to Cir X-1, assuming they are associated) to 5.5 kpc, 
which is substantially smaller than the original sug gested 
distance to Cir X-1 of 10 kpc (|Goss fc Mebold, 19771 



were found Lu dllTyi' drmiiiiLlcally eiich Lluie IL way 
observeU (see fui eAaiiiple Kalu2,ienski eL al. 1970, TenuanL 
1988, Tsunemi et al. 1989, Shirey et al. 1996). Periodic mod- 
ulation of the X-ray flux was found at a period of 16.6 d 
( Kaluzienski et al., 1976 ). A radio counterpart was found 
(Clark et al., 1975), and was fo und to flare at the s ame pe- 
riod as the X-ray modulation (Haynes et al., 1978). These 



The discovery of Type I X-ray bursts (Tennant et al. 



1986b) suggests that the compact object is probably a 



weakly-magnetised neutron star. The close association of 
Cir X-1 with the supernova remnant suggests that the sys- 
tem may be a young (< 10^ y old) runaw ay system from a 
supernova explosion (Stewart et al., 1993). 

The optical counterpart to Cir X-1 was identifled as 



flares \<| ei ' e iniiiall^i deCecCed at peak fluA le^als of > 1 Jji 
since t mo 1070o, tho flujc of tho box. 



hao dooroaaod dra 



a highly-re ddened star with strong Hq emission (Whelan 



et al., 1977). This object was later shown to consist of three 



maticall y, and it has only oc casionally been detected above 
50 mjy ( ^tewart et al., 199l| ). 

This radio source is located 25' from the centre of the su- 



stars within a radius of l'.'5, the southernmost of which is 



p ernov a remnant 0321.9 — 0.3, an d is apparently con nected 



the true counterpart ((Moneti, 1992; Duncan et al., 1993)). 

The long orbital period and periodic X-ray acti vity sug- 
gested a high-mass system in an eccentric orbit (Murdin 



to the r em nant by a radio nobula ( Haynos ot al., 1086 ). St 



et al., 1980); however, the variability of the optical emis- 



art et al. (1993) have imaged arcmin-scale coUimated struc- 
tures within the surrounding nebula, suggesting an outflow 
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sion, the faintness of the optical counterpart, and several 
of its X-ray characteristics suggest that the companion is 
a low-mass star. The lack of spectroscopic studies in the 
optical band means that most of the fundamental orbital 
parameters of Cir X-1 have not been determined. Moreover, 
Cir X-1 shows very different properties from time to time. 
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Figure 1. Spectrum of Cir X-1, taken near apastron on 1997 June 4: the spectrum shown is the sum of five 1800 s exposures. Emission 
lines of lines of Ho A6563, He I A6678, and He I A7065 can be seen; the dip at A6870 is due to imperfect removal of the terrestrial 
atmospheric B-band. 



This make it difficult to construct a coherent picture for the 
system from observations of different wavelengths at differ- 
ent epochs. Here we present new spectroscopic observations 
of Cir X-1, and use these, together with analysis of archival 
observations of the system, to suggest a more coherent model 
for the system. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 New optical observations 

Cir X-1 was observed on 1997 June 4 using the 3.9 m Anglo- 
Australian Telescope (AAT). The mean orbital phase of the 
observation wa s 0.51 , calculated according the ephemeris of 
Stewart et al. ( 1991 ). The RGO Spectrograph was used in 
combination with the TEK Ik CCD in the 82 cm camera 
and a grating of 270 grooves mm~^ in first order, result- 
ing in a dispersion of ^ 1.08 A pixel"^ over a wavelength 
ranee 6060-7165 A. The spatial scale was 0'.'25; the spectral 
resolution, measured from the arc lines, was 5.4 A. 

A l'.'5-wide slit was used, o riente d north-south so both 
Cir X-1 and star 2 of Moneti ( [l992| ) were in the slit. The 
atmospheric seeing was about 1". Five 1800 s integrations 
were taken, interspersed with CuAr arc-lamp exposures, be- 
fore cloud prevented the acquisition of any more data. 

The bias and pixel-to-pixel gain variations were re- 
moved from each exposure using standard procedures in 
I RAF. Cosmic rays were removed using the method of Croke 
(1995) to compare adjacent frames. Because of the presence 



of the nearby confusing star (Moneti's star 2), special care 
needed to be taken to measure the flux from our object. 
At every position along the dispersion direction, we flt two 
gaussian profiles, with fixed widths (FWHM=5.6 pixels) and 
separation (6 pixels), to the sky-subtracted frames. The am- 
plitude of these gaussians was used as the estimate of the 
flirx from Cir X-1 and star 2 at each wavelength. We then 
determined the wavelength calibration using the CuAr arc 
lamp exposures. We fit a low-order polynomial to the arc line 
wavelengths as a function of pixel number: the rms scatter 
of the fits was ~ 1/4 of a pixel. A rough flux-calibration was 
performed by comparing with the spectrum of the observed 
flux standard LTT 4364, though since the night was non- 



photometric, this flux calibration should be considered only 
approximate. 

2.2 Infrared observations 

K-haxiA spectroscopy of Cir X-1 was obtained using the 
Cryogenic Array Spectrometer/Imager (CASPIR) on the 
ANU 2.3 m telescope at Siding Spring Observatory on the 
night of 1997 June 20. The K grism was used with the SBRC 
256 X 256 InSb array, giving a dispersion of 2l'.'5 Apixel"^ 
over a wavelength range of 1.94-2.49 /xm. The spatial scale 
was 0'.'5 pixel" ^ A 5" slit was used, oriented east- west: note 
that this means that Moneti's stars 2 and 3 both contributed 
light in our spectrum. The telescope was nodded by ± 12" 
along the slit to provide sky frames at the same position as 
the object. Argon lamp spectra were taken to perform wave- 
length calibration, and two nearby bright stars (BS 5699 
and BS 5712) were observed in order to remove atmospheric 
spectral features and perform flux calibration. 

Standard data reduction procedures were followed, us- 
ing the local caspir package running in iraf. Bias and dark 
frames were used to linearise all frames, the sky background 
was subtracted from the object frames, and pixel-to-pixel 
variations were corrected. The chip distortion was corrected 
in order to align the dispersion and spatial directions along 
rows and columns of the chip; the sky background was then 
subtracted and spectra extracted. A low-order polynomial 
was flt to the argon lines, and these calibrations applied to 
the object spectra. Flux calibration was achieved by dividing 
the observed spectra by the spectrum of a nearby mid G- 
type star and then multiplying by a model for the absolute 
flux distribution of the calibrator. Residual terrestrial at- 
mospheric features were then corrected using an early-type 
star. 



3 RESULTS 

3.1 Optical spectra 

The optical spectrum is shown in Fig. |l]. The spectrum 
clearly shows three narrow emission lines — Ha A6563, He I 
A6678, and He I A7065 — as well as a broad component to 
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Table 1. Fit to the optical spectrum. The spectrum was nor- 
mahsed by a low-order polynomial fit to the continuum, and four 
gaussians were fit to the lines, subject to the constraint that the 
velocity and width of the three narrow lines were the same. The 
table shows the resultant velocities, FWHMs, equivalent widths 
Wx of the gaussians, and the de-reddened line luminosities, as- 
suming a reddening of Av = 11 (Predehl & Schmitt 1995) and a 
distance of 5.5 kpc (Case & Bhattacharya 1998). 



Component Velocity FWHM 
(kms"-"^) (kms~-^) 



^line 

(ergs-i) 



Ha broad 

Ho narrow 
He I A6678 
He I A7065 



-310 ±40 
-378 ±12 



2210 ±74 
402 ±8 



(A) 
74 ± 1 

3.4 ±0.5 1.5x10^"' 
4.9 ±0.5 2.1x10^"' 



the Hof line which is blue-shifted with respect to the narrow 
component. We fit gaussian models to these lines, using the 
specf it package in iraf. After normalising the spectrum by 
a low-order polynomial fit to the continuum, we fit a single 
gaussian to each of the helium lines and two gaussians to the 
Ha line. The velocities and widths of the three narrow lines 
are consistent with being the same; we therefore constrained 
them to be equal. The details of the fit are shown in Table |l| 
The narrow lines show a velocity of ~ +378 kms~^ and a 
width of ~ 9.5 A, or 400 kms~^, while the broad component 
to the Ha line has a very different velocity (—310 kms~^, 
or —688 kms~^ with respect to the narrow lines) and width 
(46 A, or 2200 kms'^). 

No stellar features can be seen in the spectrum. In par- 
ticular, no absorption features can be discerned which might 
sugggest the nature of the binary companion. There is no 
significant variability apparent between the five spectra. 

Mignani, Caraveo & Bignami (1997) obtained a spec- 
trum of Cir X-1 using the HST Faint Object Spectrograph 
near periastron (phase 0.18) in June 1995. The line profile of 
the Ha line in their spectrum was very different to our 1997 
spectrum. They found the broad component of the gaus- 
sian to be centered at approximately its rest wavelength, 
while the narrow component was observed at a velocity of 
+380 kms~^. They interpreted the narrow component as 
arising in an accretion disc, with the ~ 400 kms~^ velocity 
representing the rotation velocity of the disc material; the 
absence of the blue component they interpreted as a phase- 
dependent shadowing effect. Under this model, the velocity 
of the narrow component should shift with orbital phase, 
while the broad component remains fixed. 

Our spectrum was taken at apastron, and does not show 
the predicted shift of the narrow component. In fact, the ve- 
locities we observe for the narrow component of the Ha line 
(and the helium lines, not observed in the HST spectrum) 
are identical with the HST spectrum, while the broad com- 
ponent has shifted in velocity. This would seem to indicate 
that the model proposed by Mignani et al. for the system is 
wrong. 

To investigate this, we used data from the AAT archive, 
which contains nearly all observations taken with the AAT 
since its inception. Cir X-1 has been observed using the AAT 
several times by different observers, with most observations 
occurring about 20 years ago. We obtained spectra at a total 
of five epochs which showed detectable emission: these are 
listed in Table || and shown in Fig. |[ The May 1976 data 
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Figure 2. Line profiles of Ha in the archival observations of Cir 
X-1, showing the two gaussians fit to each line and their sum. 
Each spectrum had been binned to twice the spectral resolution, 
and normalised by a polynomial fit to the continuum, so the y- 
axis shows the relative flux of the line. The last two spectra show 
the HST observation (June 95) and our new AAT observation 
(Juno 97). 
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Table 2. List of observations of Cir X-1, including the archival AAT observations, the HST observation (June 1995), and our new AAT 
observation (June 1997). Column 2 gives the instrument/detector combination where 'B&C indicates the Boiler & Chivens spectrograph, 
'RGO' indicates the RGO spectrograph with either the long (82 cm) or short (25 cm) camera; columns 3-5 show the instrumental 
resolution, total exposure time, and orbital phase of the observation. The following columns show results of the fit to the Ha line (see 
text for details), with gaussian 1 being the red-shifted (narrow) component, and gaussian 2 the blue-shifted (broad) component. The 
phase was calculated according the ephemeris of Stewart et al. (1991), and the fits are shown in Fig. 2. 

gaussian 1 gaussian 2 

UT Date Instrument AA tcxp Phase velocity FWHM Wx velocity FWHM Wx 

(A) (s) (kms-l) (kms-i) (A) (kms"!) (kms"!) (A) 



1976 May 21 


B & C/IDS 


15 


3840 


0.536 


370 ± 40 


1050 ±40 


380 ± 40 


200 ± 60 


2760 ± 340 


200 ± 30 


1977 Feb 5 


RG025/IPCS 


0.8 


800 


0.234 


270 ± 20 


890 ± 70 


214 ±13 


-830 ±120 


1200 ±160 


57 ±10 


1977 Sep 4 


RG025/IPCS 


4.3 


2000 


0.947 


490 ± 50 


640 ± 50 


245 ± 25 


160 ±120 


2080 ± 190 


180 ± 25 


1978 Apr 29-30 


RG082/IPCS 


0.45 


10900 


0.315 


410 ±40 


720 ± 55 


52 ±6 


-190 ±70 


630 ± 80 


21±6 


1978 Aug 11 


RG025/IPCS 


1.3 


2000 


0.522 


220 ± 200 


700 ± 600 


180 ± 130 


-600 ± 2000 


1030 ±120 


30 ± 100 


1995 Jun 1 


HST/FOS 


2.2 


5160 


0.183 


330 ±10 


355 ±15 


30 ±2 


-50 ±20 


1145 ±25 


75 ±2 


1997 Jun 4 


RG082/TEK 


5.4 


9000 


0.510 


380 ±12 


400 ± 10 


74 ±2 


-300 ± 50 


2210 ±100 


50 ±2 



were published by Whelan et al. ( 1977 ) ; the rest of the data 
are unpublished. We obtained the original data from the 
AAT archive and reduced them using standard techniques 
and the FIGARO data reduction package. Note that for all 
these observations, light from stars 2 and 3 was presumably 
in the slit. Because of this, no attempt has been made to flux 
calibrate the spectra; instead, each spectrum was binned to a 
spectral sampling of two data points per resolution element, 
and the resulting spectrum was normalised by a low-order 
polynomial flt to the spectrum. 

We performed the same double gaussian fit to these 
spectra as described in Table ^ and the results are shown 
in Table |^. Several points are immediately clear; despite 
the wide range of orbital phases, the broad component 
is never observed redward of the narrow component^ 
The (weighted) mean velocity of the narrow component is 
+350 kms~^, that of the broad component — lOOkms"^. 
The scatter of the velocities of the broad component is much 
greater than that of the narrow component. No strong trends 
can be seen with either date of observation or orbital phase, 
with the exception of the equivalent width of the narrow 
component, which has been generally decreasing since 1976 
(except for the extremely low value in April 1978). 



3.2 Infrared spectra 

Infrared spectra were obtained one orbit after the optical 
spectra. The mean spectrum is shown in Fig. |^. Bright emis- 
sion lines of Br 7 and He I A2.058 were visible in the spec- 
trum. These two lines are generally observed in the K-band 



spectra of low- and high- mass X-ray binaries ( ( Bandyopad- 



hyay et al., 1997 ; Bandyopadhyay et al. 



1998 



Clark et al. 



199!:)) and may arise in accretion discs or larger emitting 



regions, such as dense winds from the mass donors and/or 
ejecta from the binaries. There also seems to be a hint of 
asymmetry in the blue wing of the infrared lines, so we again 
fit gaussians to the line profile. We modelled each emission 
line as the sum of two gaussians, constraining the widths and 
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Figure 3. Infrared spectrum of Cir X-1, taken near apastron on 
1997 June 20. 



Table 3. Fit to the infrared spectrum. The spectrum was nor- 
malised by a low-order polynomial fit to the continuum, and four 
gaussians were fit to the lines, subject to the constraint that 
the velocity and widths of the blue (B) and red (R) components 
were the same for both species. The final column shows the de- 
reddened integrated line luminosities, assuming a reddening of 
Ay = 11 and a distance of 5.5 kpc. 



Component 



Hi Br7 (R) 
He I (R) 
Hi Br7 (B) 
He: (B) 



Velocity 
( kms-l) 



FWHM 

( kms^-^) 



Wx 
(lira) 



(ergs-1) 



450 ±40 1170 ±90 



^1870 ±570 1020 ±440 



88 ±7 4.3x10^4 
96 ±7 3.1x103-* 
12 ±6 
11 ±5 



velocities of the blue and red components to be the same for 
both lines. The results are shown in Table ^ 

Both lines are well fit by one gaussian at +450 kms~*, 
a velocity very similar to the narrow red-shifted optical 
component, plus a blue-shifted component at a velocity of 
~ -2000 kms~\ 



t Except for the Apr 78 spectrum, where the width of the red- 
shifted component is marginally smaller than the width of the 
blue-shifted component; however, the two widths are consistent 
with being identical. 



4 AN ECCENTRIC LOW-MASS BINARY 
MODEL FOR CIR X-1 

Initially, it was suggested that Cir X-1 is a high mass bi- 
nary consisting of a compact star (either a neutron star or 
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Figure 4. Model for Cir X-1. During periastron (phase 0) the 
neutron star passes its closest to the companion (mass ~ 3- 
5 Mq), and the disc is disrupted, resulting in large X-ray fluc- 
tuations. After periastron passage, an accretion disc begins to 
form; mass accretion steadies, until after apastron (phase 0.5) 
steady accretion via a disc takes place. Tidal interactions begin 
to disrupt the disc again as periastron approaches. 



black h ole) and an OB supergian t companion star ((Whelan 
et al., :|977|; [Murdin et al., 198(]| )). Its X-ray properties im- 



ply a very eccentric binary orbit. The 16.6-day modulation 
in the X-ray luminosity is due to orbital variations in the 
mass accretion of the compact star. The X-ray and radio 
bursts occur when the compact star encounters the dens e 
stellar wind from the supergiant ( Fransson fc Fabian, 1980 ) , 



or when 1 :idal mass tran sfer is induced during the periastron 
passage ( [Haynes, 1987 ). The Type I X-ray b ursts that were 



discovered in a brief epoch in the mid-1980s (Tennant et al. 
1986a^ndicated that th e compact star is a ne utron star. Re- 



cent observations (e.g. ( ^tewart et al., 1991 ; Gllass, 1994)) 



indicate that the companion star is unlikely to be a massive 
supergiant star. No Type I bursts have yet been reported in 
RXTE observations. 

Here, we propose a low-mass binary model, where the 
system consists of a neutron star orbiting around a subgiant 
companion star of about 3 to 5 Mq (Fig. ^; see also Ten- 
nant & Wu 1998). The orbital eccentricity is ~ 0.7-0.9. Dur- 
ing the periastron passage, the companion star overfills its 
Roche-lobe, causing a transfer of mass at a super-Eddington 
rate onto the neutron star. As Cir X-1 is a X-ray burster, 
the magnetic field of the neutron star is relatively weak, and 
so the accretion flow is probably quasi-spherical during the 
periastron passage. Because of the large radiative pressure of 
emission from the neutron star, there must also be a strong 
(anisotropic) matter outflow. The inflow/outflow geometry 
may be related to the larger scale jets observed from the 
system, with symmetrical coUimated outflows along some 
preferred axis. 

After the periastron passage, the companion star is de- 
tached from its critical Roche-surface, and mass transfer 
ceases. Accretion continues as the neutron star captures the 
residual matter in its Roche-lobe. An accretion disc is gradu- 
ally formed. The disc is geometrically thick, because the ac- 



cretion rate is sufficiently high and is close to the Eddington 
limit. The accretion becomes more steady when disc accre- 
tion takes over from the less stable super-Eddington quasi- 
spherical accretion. Accretion via a disc continues to operate 
till the next periastron passage, where the accretion disc is 
disrupted by tidal interaction. When the periastron passage 
proceeds, the companion star overfills its Roche-lobe again, 
leading to another cycle of super-Eddington mass transfer, 
formation of the accretion disc, steady accretion via the disc 
and disruption of the accretion disc. 

The periodic X-ray light curve of Cir X-1 has changed 
dramatically since its discovery in the 70s (e.g. ( Kaluzicnski 



et al., 1976 



Tsunemi et al., 1989| ; ^hirey et al., 199^ )). This 



indicates that the accretion mode has varied substantially 
in the last 25 years. Despite the fact that the Ha lines in 
the optical spectra obtained in the 7Qs and 9Qs show certain 
similarities, such as an asymmetric profile with a broad blue 
wing, it is inappropriate to assume that the optical spec- 
tra are the same at the same orbital phases in the 70s and 
the 90s. The FWHM of the red (narrow) component of the 
Ha line seems to decrease from « 700 kms~^ in the 7Qs to 
~ 380 kms~^ in the 9Qs. The central wavelength, however, 
does not show significant velocity variation, and its average 
value is 350 ± 80 kms~^. 

Our interpretation is that the red (narrow) component 
is emitted from irradiatively-heated matter which is rela- 
tively stationary with respect to the centre of mass of the 
system and has a low velocity dispersion. One of the natu- 
ral choices is the heated surface of the companion star, on 
which the projected velocity of the matter is much less than 
the flow velocity near the neutron star.^ 

The velocity of the sharp red component of the H I Br7 
and the He I A2.058 ^m lines is 450 ± 40 kms^^, which is 
consistent with that for the optical lines. Thus, the lines 
must also be emitted from a region relatively stationary to 
the centre-of-mass of the binary. The width of the IR lines, 
however, is much larger than the width of the optical lines. 
It is therefore unlikely that the emission site is the surface 
of the companion. One possibility is that these IR lines are 
emitted from a heated dust shell - the residue of ejecta from 
previous epochs. 

The observations accumulated since the 70s show that 
the broad component is blueward of the narrow component, 
except (possibly) in the 1978 April observation. The line 
widths were generally above 1000 kms~^, and it was even 
larger than 2000 kms~^ in the 1976 and 1997 observations. 
Although the line centre velocities appear to show large vari- 
ations (from —600 kms~^ to 200 kms~^), they show no ob- 
vious correlation with orbital phase. 

We suspect that the broad component is related to the 



t Note that, if the narrow emission is arising from the heated face 
of the secondary, the equivalent width should vary in a smooth 
fashion through the binary orbit, with a maximum at phase 0.5; 
such a pattern is not seen in the data in Table jij However, we 
argue that the mass transfer rate has varied so much over the 
21 y spanned by the observations that it is misleading to try 
to compare the equivalent widths. Indeed, if we consider three 
groups of spectra near to each other in time (1976—1977, 1978, 
and 1995-1997) there is a (slight) tendency for the equivalent 
width to be highest at phase 0.5. Such a relation certainly needs 
to be confirmed. 
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Figure 5. Line profile of emission from an expanding envelope, 
with parameters as described in the text. Positive x corresponds 
to the blue-shifted frequency; negative x corresponds to the red- 
shifted frequency. The filled circles are the simulated data. The 
solid line is a fitted model with a constant and two gaussians. 
The two gaussians are represented by the dotted lines. The ratio 
of the widths of the broad and narrow gaussians is 5.5, and the 
ratio of the width of the broad gaussian to the separation of the 
two gaussian peaks is 3.22. 



high velocity fiow from regions near the neutron star. X- 
ray observations have shown that the rate of mass accretion 
onto the neutron star is above the Eddington hmit (e.g. In- 
oue 1989). It is therefore possible that an outflow is driven 
by the Eddington luminosity. As we only see the blueward 
component, the lines are probably emitted from the opti- 
cally thick outflowing matter, of which the red-shifted com- 
ponents are heavily absorbed or obscured. The changes in 
the line centre velocities and widths deduced from the fits 
are probably the result of the changes in optical depth over 
the different observational epochs. 

The quality of our data precludes us performing a fit of 
this model to the line profile; rather, we present it as a pos- 
sible explanation for the observed asymmetric line profile. 
We have constructed a simple toy model to demonstrate the 
asymmetric line profile that might be observed from such 
an outfiow. In our model the broad line emission region is 
a spherically symmetric extended envelope, with outer and 
inner radii i?out and _Rin respectively. The envelope encloses 
an opaque sphere with radius equal to i?in. We have con- 
sidered various functional forms for the velocity, emissivity 
and absorption coefficient distribution, and have found that 
models with a power-law velocity profile generally can pro- 
duce a broad line profile similar to the observations (see 
Appendix) . 

In Fig. 1^, we show a simulated a blue-shifted broad 
line from an expanding envelope with outer and inner radii. 
Rout = 1 and Rin — 0.3 respectively. The radial expanding 
speed V has a power-law profile with an index of 4-1, and 



it is normalised such that its value is 1 at i?out. The nor- 



malised local Doppler width of the line is 0.1. The emissiv- 
ity is uniform within the emitting region, with r]c = rji = 1. 
The continuum and line absorption coefficients are uniform, 
and their ratio Xi/Xc ~ 0.35. The total effective absorption 
optical depth r is such that (xc + xO^out = 1- 



As the narrow line is emitted from the heated compan- 
ion star, we simply assume it is a gaussian with a width of 
0.123 and a norm of 4.0. This will give an approximate ratio 
of 2/3 for the equivalent widths of the simulated broad and 
narrow line components. 

As shown, the resultant line profiles resembles the Ha 
line observed in our June 1997 data. If we fit the simu- 
lated line with two gaussians and a constant, the relative 
widths and shifts of the two gaussian components are in 
good agreement with the fitted parameters of the June 1997 
data (Fig. I). 

The blue components of the Hi Br7 and the He I 
A2.058 /im lines have a velocity shift of —1870 ± 570 kms~^ 
and a width of 1020 ± 440 kms~^. If we interpret the blue 
component of the IR lines as being emitted from the vicinity 
of the neutron star, the difference in the velocities between 
the optical and IR lines is due to optical depth effects and/or 
the uncertainties in determining the central velocities of the 
lines with the spect ral re solution of the IR spectrum. 

Murdin et al. ( 198C| ) assumed a 1 Mq compact star. 



and by comparing the radius of various classes of stars with 
the separation of the two component star at periastron, they 
estimated the mass of the companion star and the orbital 
eccentricity for which tidally induced mass transfer can oc- 
cur during the periastron passage. In our model, the mass 
transfer is driven by the overflow of the mass-donor star's 
Roche-lobe instead of the stellar wind from the star, and 
it therefore requires the companion star to fill its Roche- 
lobe at periastron, i.e. the star's radius R2 to be equal to 
or larger than its 'instantaneous' critical Roche-lobe sur- 
face. It is beyond the scope of this paper to derive an exact 
phase-dependent effective Roche-lobe surface. However, as 
a first approximation, if we simply take the expression of 
the Roche-lobe radius for circular systems (see Kopal 1959) 
and replace the orbital separation term by the separation be- 
tween the two stars of Cir X-1 at periastron, then we can ob- 
tain an effective Roche- lobe radius at periastron (phase 0): 



R. 



;a(l-e)[0.38 + 0.21og5] 



(1) 



where the mass ratio q = M2/M1 and a is the semi-major 
axis. Using this approximate expression for R-^^, we can con- 
strain the appropriate parameter space (e; Mi, M2) where 
Roche-lobe overflow occurs during the periastron passage for 
various types of mass-donor stars. 

In Fig. ^ we show the mass-radius relation for super- 
giant, giant and main-sequence stars, together with R-^ for 
eccentric systems with a 16.6 day period and a 1.4 Mq neu- 
tron star primary. If the companion star is a main sequence 
star, its mass is about 10 Mq and the orbital eccentricity 
e > 0.9 in order to flU its Roche-lobe at periastron. A sub- 
giant star with M2 ~ 3-5 Mq can easily fill its Roche-lobe 
for an orbital eccentricity of about 0.7-0.9. 

The argument for a low mass system is generally consis- 
tent with the combined results of our spectroscopic observa- 
tions and the study on the effects of sup ernova kicks on neu- 
tron star binaries. A simulation study (Brandt & Podsiad- 



lowski, 1995) showed that the eccentricity of HMXBs (with 
M2 ~ 10 Mq) and LMXBs (with M2 ~ 1-5 Mq) resulting 
from supernova explosions with isotropic kicks will be about 
0.5 (o" « 0.2) and 0.8 (cr « 0.05) respectively if the final or- 
bital period is about 16 days. The corresponding centre-of- 
mass velocities are 50 (a ~ 12) and 150 {a ~ 100) kms~^. 
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Figure 6. The radius of the Roche-lobe at periastron as a func- 
tion of the mass of the companion star, calculated from an ec- 
centric system with a 16.6 d period and a 1.4 Mq neutron star 
primary. For Roche-lobe overflow to occur, the size of the sec- 
ondary must exceed or equal the volume enclosed by the instan- 
taneous critical Roche-lobe surface. We also show the size of the 
star for various masses at the zero-age main-sequence (ZAMS), 
terminated main-sequence (TAMS), giant and supergiant stages 
(indicated with filled triangles, open triangles, open squares and 
open circle respectively). The radii of the stars are derived from 
the evolutionary models of Maeder & Meynet (1988), with the 
ZAMS, TAMS and giant stages corresponding to Point 1, 7 and 
14 in the evolutionary tracks. The mass-radius relation of the 
super-giant stars, which is not as well-defined, is adopted from 
Lang (1991). 



If we accept the interpretation that the narrow line compo- 
nents originate from the heated companion star, the line-of- 
sight component of the centre-of-mass velocity of the sys- 
tem is about 350 kms~^, or 430 kms~^after applying a cor- 
rection for Galactic rotation (using the rotation curve of 
Clemens 1985). This correction was made assuming that the 
position of G321.9— 0.3 was the birthplace of Cir X-1, but 
this makes little difference to the resulting correction. By 
combining with the transverse velocity (~ 300-400 kms~^) 
inferred from the age and the distan ce of the supernova rem- 
nant G321. 9-0.3 ( [Clark etal., 1975| ) we obtain an estimated 
centre-of-mass velocity of ~ 540 kms~^. This velocity could 
be substantially higher, as it depends on the (very uncer- 
tain) age estimate of the supernova remnant; transverse ve- 



locit ies as high as 1600 kms ^are possible (Stewart et al. 
1993|)7^ velocity of ~ 540 kms' 



is about 3(7 from the mean 
centre-of-mass velocities obt ained from the simulation study 



by Brandt & Podsiadlowski (1995). Such high velocities are 



not at all consistent with the simulations for high mass sys- 
tems: even the low velocity is about IScr above the mean 
velocity predicted for a high-mass system. 

The implied velocity for Cir X-1 makes it one of 
the fastest moving binaries known: radial velocities of ~ 
300 kms~^ have been me asured for a hand ful of other 
LMXBs, e.g. 4U 1556-605 (|Motch et al., 1989D, GX 349-H02 
( |Penninx fc Augusteijn, 199l| ), 4U 1957-1-11 ( [Margon et al. 
1975 ^iiee Johnston (1992) for a discussion. Such veloc- 
ities strengthen the identification of Cir X-1 as a low- 
mass system, as the fastest moving HMXBs have velocities 
V < 100 kms~^. If the transverse velocity is similar to the 



radial velocity, this implies a proper motion of 15 masyr"^; 
in an ongoing project to measure the proper motion of the 
radio counterpart to Cir X-1, we can measure a relative po- 
sitional accuracy of around 40 mas, so we should measure 
the proper motion in a couple of years. 



An alternative simulation by Shirey (1996) also con- 
cluded that Cir X-1 is a low-mass system. However, his sim- 
ulation requires the system to be old (~ 10^ y), which seems 
difficult to reconcile with the observations that Cir X-1 is as- 
sociated with the supernova remnant G321.9— 0.3. 

The RXTE observations show that the current X-ray lu- 
minosity is quite steady in the orbital phases around 0.5-1.0. 
The cycle-to-cycle variations are insignificant in comparison 
with the those at the earlier orbital phases. We interpret 
the steady accretion as due to the presence of an accretion 
disc. The erratic cycle-to-cycle luminosity variations are the 
consequence of the absence of the disc acting as a buffer; cf. 
the rms fluctuations in the RXTE data at phases 0.15 and 
0.5 (Shirey et al., 1996). Moreover, we do not see double- 
horn features in the Ha lin e, such as those often o bserved in 
cataclysmic variabl es (e.g. (Horne fc Marsh, 1986 )) or black 
hole binaries (e.g. ( 3oria et al., 1998| )). This, together with 
the near-Eddington luminosity of the X-rays during the disc- 
accretion phase, implies that the disc is non-Keplerian and 
it is both geometrically and optically thick. As the matter 
is not tightly bound by the gravity of the neutron star, the 
large variation in the tidal force when the neutron star ap- 
proaches periastron could easily disrupt the disc, causing a 
short temporal decrease in the accretion luminosity of the 
system. 

What constraints can we put on the age of the system? 
In our model, mass transfer at periastron does not begin un- 
til the star has evolved from the main-sequence towards the 
giant branch. For a 3-5 Mq star, the subgiant stage may 
last ~ 5 X lO' y. If Cir X-1 is associated with the nearby 
supernova remnant G321— 0.3, then the implied age of the 
system as an eccentric neutron star binary is < 10^ y. The 
Eddington limit for the rate of mass transfer onto a neutron 
star is ~ 3 X 10* Mq y~^, giving a mass transfer time scale 
of 3 X 10^ y. Since the companion star in Cir X-1 can lose 
mass at a rate greater than this as an average in each or- 
bital cycle, the life-span of the system as an X-ray binary is 
correspondingly shorter, but still comfortably above either 
the evolution time scale of the companion star or the age of 
the supernova remnant. 



5 CONCLUSIONS 

We have detected strong optical and infrared emission lines 
in the spectrum of Cir X-1; the Hq and infrared lines 
are asymmetric, with a narrow component at a velocity 
of ~ -1-350 kms~^ and a broader, blue-shifted component. 
Archival optical observations show that an asymmetric Ha 
emission line has been in evidence for the past twenty years, 
although the shape of the line has changed significantly. The 
narrow component is always seen redwards of the broad com- 
ponent, at a velocity around 200-400 km s~^ ; the broad com- 
ponent has a much larger scatter of velocities. These prop- 
erties a re no t consistent with the interpretation of Mignani 
et al. (1997), who suggested that the narrow component 
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arises in an accretion disc, with the absence of a red com- 
ponent due to a phase-dependent shadowing effect. 

Instead, we suggest that the narrow component arises 
from the heated surface of the companion star, which is 
a subgiant of about 3-5 Mq, while the broad component 
arises in an optically thick, high velocity outflow driven by 
super-Eddington accretion onto the neutron star. During 
periastron passage, the companion star overfills its Roche- 
lobe, causing a transfer of mass at a super-Eddington rate, 
which in turn drives a strong matter outflow. After perias- 
tron, mass transfer from the companion ceases, but accre- 
tion continues at a near-Eddington rate as the neutron star 
captures the residual matter in its Roche-lobe. An accretion 
disc gradually forms, which is non-Keplerian and geomet- 
rically thick. This change between quasi-spherical and disc 
accretion causes the change between strong X-ray variability 
near phase and steady accretion between phases 0.5-1.0. 
Thus our observations are in general consistent with a low- 
mass binary model for the system. In this model, the veloc- 
ity of the narrow component reflects the space velocity of 
the binary; the implied radial velocity (-1-430 kms~^) makes 
Cir X-1 one of the fastest moving binaries known. 
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APPENDIX A: LINE PROFILE SIMULATIONS 

The transfer equation that determines the intensity of the 
radiation is 



cdt 



+ (n- V) + x(«^,r) 



I{u,r) = ??(i',r). 



(Al) 



where I{i',r) is the intensity, h is the unit vector in the 
direction of the ray propagation, ri{i/, r) is the emissivity 
and x(i', is the extinction coefficient. We assume a (quasi- 
)stationary state, such that the term d/dt can be set to zero. 
In a general situation, the extinction coefficient consists of an 
absorption term and a scattering term. Therefore, one needs 
to solve an integro-differential equation to determine the 
radiation field and then the emerging intensity. However, if 
we can neglect scatterings, the situation is greatly simplified 
and a simple formal solution can be obtained. 

In our model, we assume that the radiation is emitted 
from an envelope within which matter is moving with veloc- 
ity v{r). Because the emitters are in relative motion to the 
observer, there are velocity shifts in the frequencies of the 
radiation they emit. There is also a velocity spread among 



-2 



Figure Al. Simulated broad line components from expanding 
envelopes for a power-law radial velocity profile with various 
power-law indices. From top to bottom, Xl/Xc = 0-1) 1 10. 
The other parameters are: 77; = 1, ?7c = 1, -Rout = 1, ^?in = 0.1, 
5 = 0.1 and (x; + Xc)Ro\it = 1. The velocity normalisations corre- 
sponding to the three cases with the power-law indices -|-1 (solid 
line), (dotted line) and —1 (dashed line) are ^(-Rout) = 1, 1 
and 0.45 respectively. 



the emitters, and so the frequency-dependence of the local 
line absorption and emission coefficients are modified. 

For mathematical convenience we consider a dimension- 
less frequency x{= c{iy — 1^0) /voV*), and a dimensionless ve- 
locity V(f){= v{r)/Vf,), where i^o is the frequency in the 
rest frame, w* is a velocity normalisation and c is the speed 
of light. The emissivity and the absorption coefficients are 
decomposed into a continuum and a line component. As a 
first approximation, we assume that the continuum compo- 
nent does not vary significantly with frequency near the line 
frequencies, and hence we have 



r]{x,r) = rjc{r) +r]i{r)(t){x,f,V) ; 
X{x,r) = Xc{r) +Xi{r)4'{x,r,V) . 



(A2) 
(A3) 



In the above expression, 
by 

1 



{x, f, V) is the line profile, given 



ix,r,V) = 



(5(r)i/7r 



exp{- 



(A4) 



where 5{r) is the normalised local Doppler width. 

As the rays that reach a distance observer are parallel, 
we can assign each of them an impact parameter p with 
respect to the centre of the emitting envelope. The flux at a 
frequency x seen by the observer is the sum of all rays, i.e. 



F{x) — dip dppl*{x,p) 



(A5) 



where I*{x,p) is the intensity of the ray emerging from non- 
obscured part of the emitting envelope 5R, and J dip J dp p 
is the projection area of the expanding envelope on the sky 
plane. If we neglect scatterings, I*{x,p) is the formal solu- 
tion to the radiative transfer equation, and it is 



r{x,p) 



Io{x,p) exp{~T{x,p,Zmin{p))} 



(A6) 
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k(p) 



a{p) 



dz [i7c + ??!0] exp{-r(a;,p,2:(p))} 



(e.g. dMihalas, 1978| [Kudritzki et al., 1988[ )), where dz is 
a path element along the ray propagation, and Io{x,p) is 
the intensity of the background radiation. The optical depth 

t{x,p,z{p)) is 

(•Zmax{p) 

t{x,p,z{p)) = / dz[xa + Xi'f\ ■ (A7) 

Jz(p) 

If there is no strong background radiation behind the emit- 
ting envelope (i.e. Io{x,p) ~ 0), then the observed flux is 
simply 

d(p / dpp / dz[ric + Tjicf^ 

xexp{-T{x,p,z(p))}. (A8) 

We consider a hybrid method to calculate the line flux 
F{x). The optical depth t{x,p,z{p)) is determined by a di- 
rection integration of the absorption coefficient along the ray 
propagation hz (for fixed p). Instead of carrying out a direct 
integration in the 3-D space (v?,p, z), we use a Monte-Carlo 
method to determine F{x). 

In the Monte-Carlo calculation, we consider each emit- 
ter as a discrete source which emits a packet of line and con- 
tinuum radiation spanning a certain frequency range. The 
packet is absorbed along the path of propagation subject to 
the optical depth. The strength of the emission from each 
emitter (i.e. 77c and rji) is weighted according to a distribu- 
tion profile in the envelope. We consider 50000 emitters in 
each simulation. The observed flux is the sum of the contri- 
bution of all emitters in the non-obscured region 3?. 

This quick-fix method allows us to handle complicated 
geometry due to obscuration easily. It also enables us to 
survey various functional forms for the distribution profiles 
of the velocity, emissivity and absorption coefficients and 
obtain qualitative features of the line emission with sensible 
computational time. 

In Fig. Al we show examples of line profiles for three 
cases in which the radial velocity distribution is a power-law 
function (with indices —1, and -1-1 respectively). In these 
cases, if the line absorption is weaker than the continuum 
absorption, i.e. xi < Xc, an asymmetric blue-shifted profile 
can be obtained. If xi ~ Xc, the line is symmetric with no 
obvious shift of the line centre frequency. If Xi > Xc, a P- 
Cygni profile will be obtained. 
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